Theory of tunneling spectroscopy in the Larkin-Ovchinnikov state 



O 
O 



(N 
O 

o 

Oh' 
c3 



^3 

o 
o 



in 

O 

o 

O 

-i— > 
c3 



i 

c 

o 
o 



X 



Y. Tanaka 1 ' 2 , Y. Asano 3 , M. Ichioka 4 and S. Kashiwaya 5 
1 Department of Applied Physics, 
Nagoya University, Nagoya, 464-8603, Japan 
2 CREST, Japan Science and Technology Corporation (JST) Nagoya, 

464-8603, Japan 
3 Department of Applied Physics, 
Hokkaido University, Sapporo 060-8628, Japan 
4 Department of Physics, Okayama University, 
Okayama 700-8530, Japan 
NeRI of National Institute of Advanced Industrial Science and Technology (AIST), 

Tsukuba, 305-8568, Japan 



A theory of the tunneling spectroscopy of normal metal / Larkin-Ovchinnikov (LO) state is 
presented by fully taking account of the periodic modulation of pair potentials in real space. The 
resulting tunneling spectra have characteristic line shapes with several maxima and minima reflecting 
minigap structures due to the periodic pair potentials. These features are qualitatively different from 
tunneling spectra of Fulde-Ferrell state and those of uniform s- and d-wave superconductors. On 
the basis of calculated results, we propose a method to identify the LO state in actual experiments. 



Fulde and Ferrell and Larkin and Ovchinnikov 
proposed superconducting states in which spin-singlet 
superconducting pair potentials are periodically modu- 
lated in real space under high magnetic fields. The Zee- 
man spin splitting results in a total momentum 2q of 
Cooper pairs. Fulde-Ferrell (FF) discussed that pair po- 
tential becomes Aexp(iq • r). Larkin and Ovchinnikov 
(LO) proposed independently an alternative scenario in 
which the order parameter is real, but varies periodi- 
cally in real space like Acos(q • f) 0- The two states 
are collectively called Fulde-Ferrell-Larkin-Ovchinnikov 
(FFLO) state. Only a few attempts have been made 
at FFLO state in the last century. However FFLO 
state suddenly becomes a hot topic because evidence of 
FFLO state has been reported in a number of experi- 
ments 0, 0, IE an d theories 0, . These studies 
suggest realization of FFLO state in a heavy fermionic 
compound CeCoIn.5. According to the phase diagram of 
CeCoIn 5 , LO state is considered to be more stable than 
FF state. To understand basic properties of LO state, it 
is important to study detailed structures of quasiparticle 
energy spectra. In fact, experimental research in this di- 
rection becomes accessible now 0, ^| . Thus a theory 
of the tunneling spectroscopy of LO state is desired to 
interpret the tunnel spectra in experiments. 

In the presence of the spatial modulation of pair po- 
tentials 0, Q], local density of states (LDOS) is ex- 
pected to be very different from that in the conventional 
[/-shaped subgap structure in uniform superconductors. 
Up to now, however, there is no theory of tunneling spec- 
troscopy in which effects of the periodic sign change of 
pair potential in real space on spectra are fully taken into 
account. In addition, wc also have to consider effects of 
sign change of pair potential in fc-space because CeCoIns 
is a strongly correlated material and its promising pair 



potential is d-wave symmetry. In such a pairing symme- 
try, charge transport of the junctions are governed by a 
mid gap Andreev resonant state (MARS) formed at the 
junction interface [IE EE E3, El EH- Thus we must 
take into account the two kinds of sign change in pair 
potentials, i.e., the sign change in real space and that in 
fc-space, on an equal footing. 

In this paper, we present a theory of tunneling spectra 
in normal metal / insulator / superconductor in LO state 
(N/I/LO) for the first time. We choose that q is parallel 
to the normal of the junction interface. In such junctions, 
tunneling spectra are expected to have fine structures re- 
flecting the miniband due to the spatial modulation of 
pair potentials. We consider two types of junction: (1) 
node contact junctions where the amplitude of LO pair 
potentials vanishes at the junction interface and (2) max- 
imum contact junctions where the amplitude of pair po- 
tentials takes its maximum at the junction interface. In 
addition to tunneling spectra, we calculate LDOS of bulk 
LO state at two local points: one is a place at which pair 
potential takes its maximum (maximum point) and the 
other is a place at which pair potential has a node (nodal 
point). Calculated results show complex structures in 
tunneling spectra and LDOS depending on types of junc- 
tion and pairing symmetries. These are a consequence of 
characteristic quasiparticle excitation in LO state such as 
miniband structures, Zeeman spin splitting and MARS 
at the interface. By comparing tunneling spectra and 
LDOS, we discuss a relation between quasiparticle states 
originating from the sign change of pair potential in k- 
space and those originating from the sign change in real 
space. The calculated results give us useful information 
to identify the pairing symmetry of LO state in CeCoIns . 

Let us consider an N/I/LO junctions in two-dimension 
as shown in Fig. 1(a), where N and LO correspond to 
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regions for x < and x > 0, respectively. We assume a 
flat interface in the y direction and isotropic Fermi sur- 
face. Magnetic field H is applied in the x direction so 
that effects of magnetic fields on orbital part of wave 
function can be neglected. The insulating barrier is ex- 
pressed by the delta-function model and its potential is 
given by Hf,5(x). The effective mass m and Fermi wave 
number kp are chosen to be common in N and LO. Here, 
we neglect the spin-orbit coupling in LO. The tunneling 
conductance in the junctions 175 at zero temperature can 
be expressed by the summation of up and down compo- 
nents with 
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with V being a bias voltage [1?J. The transparency at 
the interface is given by a n (9) =4 cos 2 9 /(4 cos 2 9 + Z 2 ), 
where Z = 2mHi,/hkp x , 9 is an injection angle of a quasi- 
particle at the interface measured from the x direction 
and Ufx is the x component of the Fermi wave vector. In 
the above, r±j and T±i obey the following equations 

-it> Fa x r ±T =A(0±,a;)(l + 1%) - 2r ±T (e + ^), (4) 

-iv Fx d x T ±l =A(0±,x)(l + r 2 ±i ) - 2r ±i (e - ^), (5) 

where 9 + = 9, 9_ = n — 9, e denotes a quasiparticle en- 
ergy measured from the Fermi energy, ujl is the spin Lar- 
mor frequency, and vp x is the x component of the Fermi 
velocity, respectively. The spatial dependence of pair 
potential A(6,x) is given by A(6,x) = A(x)Q(x)f(9), 
where Q(x) is a step function. The form factor f(9) 
depends on pairing symmetries and orientations of crys- 
talline axis in CeCoIns: f(ff) = 1 for s-wave symmetry, 
f{9) = cos26> for d a .2_ J ,2-wave symmetry, and f(9) = 
sin 29 for d xy -wetve symmetry. We calculate r±jm i n 
Eq. J3J from 2x2 quasiclassical Green's functions based 
on the evolution operator method |2l|. We choose the 
value u>l as wl = OSttTq with Tq being the transition 
temperature. The resulting effective magnitude of the 
Zeeman splitting is about 0.7Ao with Ao ~ 0.56irTc- In 
the present paper, we describe the spatial dependence of 
LO pair potentials by A(x) = AoCOs(Qx + Qx), where 
Q = 2it/L with L being the period of oscillations in 
pair potentials. The period L is measured in units of 
Lq = 2tt /Qq with Qo = The barrier parameter Z 

at the interface is chosen to be Z = 5 throughout this pa- 
per, and resulting transparency of junction is about 0.1. 
In what follows, we discuss normalized tunneling con- 
ductance ot = o's(eV)/aN, where &n — J d9 cos9a n is 
the normal conductance of junctions. In addition to <jt, 
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FIG. 1: A normal metal / insulator / superconductor in LO 
state junction is shown in (a). In a node contact junction and 
a maximum contact junction, Qi is set to be —tt/2 and 0, 
respectively. A schematic illustration of scanning tunneling 
spectroscopy is in (b), where a cone represents a scanning 
probe for local density of states. 



we also calculate pt which is LDOS normalized by its 
value in the normal state. In experiments, LDOS can be 
measured in scanning tunneling spectroscopy as shown 
in Fig. 1(b). 

First we focus on the s-wave symmetry. In Fig. 2(a) 
and (c), (Xt of the node contact junction and that of the 
maximum contact junction are shown, where we choose 
period of oscillations to be L = Lq for dotted line a and 
L = 10L for solid line b. In Figs. 2(b) and (d), px at 
a nodal point and that at a maximum point of bulk su- 
perconductors are shown. The LDOS at the nodal and 
maximum points of the pair potential are significantly 
different from each other [20( as shown in Figs. 2(b) and 
(d). In both types of junctions, ctt and px have fine 
structures reflecting the miniband structures due to spa- 
tial oscillations of pair potentials. When we compare 
curve a with curve b, fine structures are more remark- 
able in longer period of oscillations. In the maximum 
contact junctions, line shapes of ut and those of pp are 
very similar to each other. This feature becomes clearer 
for longer period of oscillations. On the other hand in 
the node contact junctions, line shapes of <7t are very 
different from those of px- Here we note that results 
for (4,2 _„2 -wave symmetry are qualitatively the same as 
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FIG. 2: Calculated results for s-wave N/I/LO junctions. Nor- 
malized tunneling conductance ctt is plotted as a function of 
bias voltages for the node contact junction in (a) and for the 
maximum contact junction in (c). Local density of states pr 
at a nodal point and that at a maximum point of bulk LO 
state are shown in (b) and (c), respectively. In all panels, a: 
L — Lo and b: L — IOLq. 



those in the s-wave junctions in Fig. 2. 

Secondly we focus on e^-wave N/I/LO junctions. In 
this case, MARS is formed at the interface due to the sign 
change of the pair potentials in fc-space. When spatial 
modulation of pair potentials is absent, ctt has a splitting 
peaks around the zero bias. The splitting width corre- 
sponds to the Zeeman splitting energy. We show calcu- 
lated results in Figs. 3(a)-(d). The line shapes of the ot 
in the maximum contact junction in (c) are very much 
different from those of px in (d). On the other hand in 
the node contact junctions, the line shapes of the ctt in 
(a) are similar to those of px in (b) . The correspondence 
between <jt and pr is tending to be clearer for the longer 
period LO state. A number of peaks and dips due to the 
miniband structures can be seen in curve b in Figs. 3(a) 
and (b). 

We summarize obtained results in Figs. 2 and 3, and 
explain the characteristic features by using schematic pic- 
tures in Fig. 4. The sign change of pair potential in k- 
space affects drastically ot and that in real space mod- 
ifies pt- In the maximum contact junctions of s-wave 
symmetry, ot well corresponds to pr because the sign 
change of pair potential at the interface is absent in both 
fc-space and real space. On the other hand in the node 
contact junctions of (icy-symmetry, <jt well corresponds 
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FIG. 3: Calculated results for d xy -wa,ve N/I/LO junctions. 
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FIG. 4: Schematic illustration of the sign change of the pair 
potential. Correspondence between the reflection process at 
the interface and the local density of states in the bulk LO 
state. 



to pt- This is because <jt and pr reflect the sign change 
in fc-space and that in real space, respectively. Thus tun- 
neling spectra are closely related to LDOS even in the 
presence of MARS at the interface. These features are 
specific to LO state. 

Thirdly we look into ot for FF state as shown in Fig. 5. 
In FF state, pair amplitudes are uniform in real space. 
Thus a preexistence formula for uniform superconduc- 
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FIG. 5: Calculated results for N/I/FF junctions. Normal- 
ized tunneling conductance ar is plotted as a function of bias 
voltages for s-wave in (a) and that for d^-wave in (c). Lo- 
cal density of states pr for s-wave and that for d xy -wme are 
shown in (b) and (c), respectively. 



tors can be applied to the N/I/FF junctions even in 
a presence of the Zeeman splitting. The calculated re- 
sults of pt do not have a spatial dependence contrary to 
those of LO state. For s-wave case, the line shape of <jt 
becomes similar to that of px for longer period FF state 
(see curve b in (a) and (b)). On the other hand, for d xy - 
wave case, the line shape of <tt has a splitting zero-bias 
conductance peak (ZBCP) due to the sign change of pair 
potentials in fc-space |22|. For almost all of the cases, 
the overall features are very different from those in LO 
junctions. Only for the case of FF state with d xy pairing, 
the line shape of ot of longer period junctions is similar 
to the corresponding case of LO state. 

Finally we propose a tunneling experiment at a suf- 
ficiently low temperature and predict observed tunnel 
conductance on the basis of the obtained results. Since 
CeCoIns is a strongly correlated electron system, real- 
ization of the (i-wave state is more plausible than that of 
s-wave [13, 13- ^ th e (HO) crystalline axis of CeCoIn 5 
is perpendicular to junction interface, which corresponds 
to d X y-~wave symmetry in this paper, a ZBCP is expected 
in tunneling spectra in the absence of magnetic field. We 
note that ZBCP can be always expected in <i-wave junc- 
tion except for an idealistic junction in which (100) crys- 
talline axis of CeCoIns is just perpendicular to junction 
interface. At the zero magnetic field, ctt has a ZBCP. 
The height of ZBCP would decrease with the increase of 



magnetic fields and the ZBCP would split into two peaks 
because of the Zeeman effect. The splitting width of 
the ZBCP increases with the increase of magnetic fields. 
When we increase magnetic fields further, CeCoIn5 un- 
dergoes transition to the LO state. After the transition, 
fine structures with many peaks and dips are expected in 
ctt in addition to the splitting ZBCP. When CeCoIn 5 un- 
dergoes the transition to the normal state at the critical 
magnetic field, the splitting ZBCP and the fine structures 
would disappear. If above change in tunneling spectra 
would be observed, it must be concrete evidence of the 
LO state. 

In conclusion, a theory of the tunneling spectroscopy 
of normal metal / Larkin-Ovchinnikov (LO) state is pre- 
sented by fully taking account of the periodic modulation 
of pair potentials in real space. The resulting normalized 
tunneling conductance has a line shape with several max- 
ima and minima reflecting minigap structures of density 
of states. These features are not expected in the Fulde- 
Ferrell (FF) state because amplitude of pair potentials is 
uniform in the FF state. Our results serve as a guide to 
identify the LO state in the experiments. 
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